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ABSTRACT 

The solar surface dynamo has become an active area of research in an attempt 
to understand the origin of a variety of magnetic field structures on the sun. 
The major modification that needs to be incorporated in the standard dy- 
namo process is the inclusion of the partial ionization of the gas in the layers 
underlying and overlaying the photosphere along with the effects associated 
with the multifluid nature of the system. This not only changes the inertia 
carrying species but also substantially modifies the temporal and spatial evo- 
lution of the magnetic induction. The energy equation also carries the import 
of these non-ideal effects. The effects such as the Hall effect and the ambipolar 
diffusion take the dynamo study beyond the realm of the ideal magnetohy- 
drodynamics. In this paper, a first principle formulation of the solar surface 
dynamo problem has been attempted. 

Key words: partially ionized plasma, surface dynamo, Hall effect, ambipolar 
diffusion . 



1 INTRODUCTION 

Vogler and Schiissler ( 2007) with their paper on solar surface dynamo have opened up a new 
area of the investigation of the solar magnetic field structures. The solar surface fields on 
relatively short spatial and time scales with varying geometries are now observed and inferred 
routinely. This necessitates the inclusion of the realistic solar atmospheric models in the 
dynamo mechanism. Vogler and Schiissler ( 2007) have studied the surface dynamo including 
the effects of stratification, compressibility, partial ionization, radiative transfer as well as 
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an open lower boundary. While this is a major effort, the magnetic induction equation they 
have used is still bereft of the effects ensuing from the multifluid interactions in a partially 
ionized plasma specifically the Hall and the ambipolar diffusion effects. Recently Krishan and 
Gangadhara (2008) initiated a study of the meanfield dynamo in a partially ionized plasma 
including the Hall and the ambipolar effects in the induction equation for closed boundary 
conditions while neglecting stratification, compressibility and radiative transfer effects. In 
this paper an attempt is made to formulate the solar surface dynamo including the Hall effect 
and the ambipolar diffusion while retaining all the effects included in Vogler and Schiissler ( 
2007). A single fluid description of the partially ionized plasma in the limit of weak ionization 
is presented in section 2. The mean-field dynamo equation in the turbulent partially ionized 
plasma is derived in section 3 and the expressions for the a and the (5 parameters of a 
meanfield dynamo are extracted. A crude estimate of the time scale of evolution of the 
magnetic field due to turbulent ambipolar diffusion arising due to the fluctuating ion-neutral 
collision frequency is given. 



2 THREE-COMPONENT MAGNETOFLUID 



We begin with the three component partially ionized plasma consisting of electrons (e), 
ions (i) of mass density pi and neutral particles (n) of mass density p n . A weakly ionized 
plasma is defined by the condition (Alfven and Falthammer 1962) that the electron-neutral 
collision frequency u cn ~ 10~~ 15 n n \j8K B T/(irm ea ) is much larger than the electron- ion colli- 
sion frequency u ci ~ 6 x lO~ 24 niAZ 2 (K B T)~ 3 / 2 . This translates into the ionization fraction 
n e /n n < 5 x 10~ n T 2 (Alfven and Falthammer 1962) where n's are the particle densities and 
T is the temperature in Kelvin. A major part of the solar photosphere (Leake & Arber 2006; 
Krishan & Varghese 2007) qualifies to be a weakly ionized plasma. The equation of motion 
of the electrons can be written as: 
~dV 



m e n e 



dt 



+ (K • V)v c 



erir, 



E + 



V x B 



= -Vp e - 

m e n e v en (V e - V n ) . 



(1) 



where the electron-ion collisions have been neglected since the ionized component is of low 
density. On neglecting the electron inertial force, the electric field E is found to be: 
V c x B Vp e m e 



E 



V~ 



- V n ) . 



(2) 
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This gives us Ohm's law. For 5 = (pi/p n ) <C 1 along with the dominance of the ion-neutral 
collisional force over the inertial force, the ion dynamics can be ignored. The ion force balance 
then becomes: 

Vi x B 



= — Vpi + erii 



E + 



- V n ) , 



(3) 



where v m is the ion-neutral collision frequency, and the ion-electron collisions have been 
neglected for the low density ionized component. Substituting for E from Eq. (2) we find 
the relative velocity between the ions and the neutrals: 
V(pi +p c ) J x B m c u en 



where 



J, 



J = en c (M - V e ) 



(4) 



(5) 



and a term proportional to (m c nu cn ) / (miU^) has been neglected. The equation of motion of 
the neutral fluid is: 
~dV n 



Pn 



dt 



+ (K-V)V n 



= -Vp n - V n iPn(V n ~ V,) - 

^nePn(K ~ V e ) + p n g + 
V.T, 



(6) 



where g is the gravitational acceleration, r is the viscous stress tensor of the neutral fluid 
defined as 



'dV r 



in ^^nj ^ 



^ (V.V n ) , 



i,j = 1,2,3. . 



Substituting for "V^ — "V^ from Eq. (4), and using v- m p x = h , ni p Q we find: 



(7) 



Pn 



OK 

dt 



+ (V n -V)K 



„ J x B 

-Vp + + p n + 

C 

V.r, . 



(8) 



where p = p n + Pi + Pc and /i is the dynamic viscosity. Observe that the neutral fluid is 
subjected to the Lorentz force as a result of the strong ion- neutral coupling due to their 
collisions. A comparison with the corresponding equation of Vogler and Schiissler ( 2007) 
shows that here the neutral mass density p n appears instead of the total mass density. This 
is valid for p n » pi in a weakly ionized plasma. 
Consider Faraday's law of induction: 



dB 

~dt 



-cV x E 



(9) 
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By substituting for the electric field from Eq. (2), we get 
dB 
~dt 



V x 



V e xB + -^V ft +^(V,-V n ) 

en P e 



(10) 



Here 77 = m e z/ en c 2 /(47re 2 n e ) is the electrical resistivity predominantly due to electron- 
neutral collisions. Using the construction 

V c x B = [V n - (V n -Vi) - (Vi- V e )} x B , (11) 

substituting for the relative velocity of the ion and the neutral fluids from Eq. (4), assuming 
p c = pi and neglecting terms of the order (m e u en ) / (miU^) Eq. (10) becomes: 

V x 



dB 



dt 



J JxB\ _ 

V n + x B + 

en e cu- m pi J 



r]V 2 B . (12) 

where S~2i C = ^ is the ion cyclotron frequency. One can easily identify the Hall term (J/en e ), 
and the ambipolar diffusion term (J x B) (Brandenburg and Zweibel 1994, Chitre & Krishan 
2001). The Hall term is much larger than the ambipolar term for large neutral particle 
densities or for u m ^> u ci . In this system the magnetic field is not frozen to any of the fluids. 
The induction equation carries the effects of stratification as well as the Vp; x B fluid drift. 
In what follows the electron and the ion densities (pressures) are assumed to be equal and 
constant. We write the induction equation in a compact form as: 



OB 



= V x [V E x B -77V x B] , 



dt 
where 

v E = v n + v H + V Am 

with 

en e 

as the Hall velocity and 
JxB 

'Am — 



(13) 



(14) 



(15) 



(16) 



could be called the ambipolar velocity. One can notice that the Hall effect and the ambipolar 
diffusion terms are not included in the induction equation given by Vogler et al. (2005) and 
used by Vogler and Schiissler (2007) to study solar surface dynamo. In order to appreciate 
the time and the spatial scales over which the additional terms would contribute to the 
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magnetic induction evolution, the induction equation is written in a dimensionless form by 
using the normalizing quantities B , L ,to,Vo,n n0 such that Vq = L /t = (Bl/Aitpn) 1 / 2 . 
The dimensionless induction equation reads: 



OB 

~dt ~ 

where 



V x [(V n - e H W x B + e A (V x B) x B) x B - r?V x B\ , 



L o \Pi, 



1/2 



(17) 



(18) 



and Xi = is the ion inertial scale and uj iv = ( 47T ^ e ) 1 ^ 2 is the ion plasma frequency. 
Thus the Hall term brings in a finite spatial scale in an otherwise ideal MHD system of any 
arbitrary scale L. The coefficient of the ambipolar term is found to be: 



£a = — £h ■ 



(19) 



clearly highlighting the predominance of the ambipolar diffusion over the Hall effect for large 
B and small v in . 

The mass conservation is described as: 

^ + v-( Pn y n ) = o, 
v-m = o, 
v • v c = o . 

The energy equation for a partially ionized plasma is found to be: 



(20) 



<9e ^ 



+ (Vh + V Al 



B' 



[B.V F )B 



-V. 



8tt 4tt 
>(V x B) x B 

47T 



+ 



[V. (V n .r) - 8^] + V. («VT) + 

p n (g-V n ) 



^J 2 + Q. 



(21) 



where Sj,- = Here Q includes the rest of the nonadiabatic effects such as the radiative 
cooling. The energy equation explicitly displays the additional terms due to the Hall and 
the ambipolar effects. The effect of the ionization potential \% nas n °t been taken explicitly 
into account in the internal energy. Instead one could continue with the methodology used 
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by Vogler et al. (2005) to consider pressure p and temperature T to be a function of the 
internal energy . 

Equations (8), (13), (20) and (21) describe the solar surface dynamo. 

3 THE ALPHA EFFECT IN COMPRESSIBLE THREE-COMPONENT 
TURBULENT MAGNETOFLUID 

The a effect in an incompressible partially ionized plasma has been discussed recently ( 
Krishan and Gangadhara 2008). Here an attempt is made to include the effect of turbulent 
density fluctuation of the neutral fluid along with the usual velocity and magnetic field fluc- 
tuations. The ambipolar diffusion term depends on the density of the neutral fluid through 
the ion-neutral collision frequency. The electron and the ion densities are kept constant and 
their fluctuations are ignored. Following the standard procedure (Krause & Radler 1980) 
the velocity Ve, the magnetic field B and the neutral fluid density p n are split into their 
average large scale parts and the fluctuating small scale parts as: 

V E = Ve + VI (22) 
B = B + B', (23) 

Pn = Pn~ + P n (24) 

such that 

n = 0, W = 0, ^ = 0. (25) 

In the mean-field dynamo the magnetic induction equation is solved for large and small 
scale fields. Here, the mass conservation equation of the neutral fluid would also be consid- 
ered. Substituting Eqs. (22-25) into the induction equation (13), we find, in the first order 
smoothing approximation, 

K = K -^ + ^- + ^-^lxB. (26) 
en c cz/ in pi cz/ in pi cpiU in z 

and the mean flow is found to be: 



_ _ J JxB J'xB' JxB'v n 
V E = V n + + — = =^ 

en e cu in pi cu in pi cpiV in A 



J 1 x Bis' 



2 

CPi Vin 



(27) 



The fluctuation in the ion-neutral collision frequency is defined as: 

v' m = Ap' n , . (28) 
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i^M" 1 . (29) 
m n \TTm in J 

where the cross-section S in ~ 5 x 10~ 15 cm 2 and m in = m ' l " 1 " is the effective mass of the 

"* mi+m n 

scattering particles. 

The turbulent electromotive force S is a function of the mean magnetic induction B and 
mean quantities formed from the fluctuations, and is expressed as: 



S = V£ x B' = aB-(3V x B , (30) 
where 



« = -y^'(Vx^) 

= a v + a H + «Am • (31) 
Here 



«v = -^K'^„ (32) 

is the measure of the average kinetic helicity of the neutral fluid in the turbulence possessing 
correlations over time r cor . Retaining only the first order contributions from the Hall and 
the ambipolar effects, the contribution to a from the Hall effect is found to be 

2r„ 



oh = ir^J' ■ K (33) 



The coupling of the charged components with the neutral fluid is clearly manifest through the 
possible correlation between the current density fluctuations and the vorticity fluctuations 
of the neutral fluid Q' n — V x V^. The ambipolar term gives rise to 

«Am = CKA ' B + CK Al • J + «A 2 J X B , (34) 

with 

«a = J^K , (35) 

3cpiZ/ in 

<*M = --^S^T^ , (36) 
3cpiZ/ in 



and 



A\ 2r„ 



<*a 2 = (— ^'xVA. (37) 
V z/ in / 3cpiZ/ in 

The contributions to a from the ambipolar diffusion with its essential nonlinear character 
manifest through its dependence on the average magnetic induction and its spatial variation. 
Note that the term okai was ignored in Krishan and Gangadhara (2008) as therein only B 
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dependent terms were retained. Now that the terms proportional to J x B appear, it becomes 
necessary to retain terms proportional to J. The compressibilty effect appears through the 
correlation between the gradient of the density fluctuation and the velocity fluctuation and 
this term brings in a quadratic dependence on the large scale field B. One also observes 
that the Hall alpha (Eq. 33) requires a component of the fluctuating current density along 
the fluctuating vorticity of the neutral fluid whereas the ambipolar effect (Eq. 35) thrives on 
the component of the fluctuating current density perpendicular to the fluctuating vorticity. 
The turbulent dissipation parameter (3 is given by 



P = ^ L V£=p v + p H + p Am (38) 
with 



Pv = ~y v n 2 (39) 

as the measure of the average turbulent kinetic energy of the neutral fluid in the turbulence 
possessing correlations over time r cor . The first order Hall contribution to (3 is 

Pu = -P^^K (40) 
3en e 

The coupling of the charged components with the neutral fluid is clearly manifest through 
the possible correlation between the current density fluctuations and the velocity fluctuations 
of the neutral fluid. The ambipolar term furnishes 

p Am = (3 A0 -B + f3 A1 -J + (3 A2 -(jxB) , (41) 
where 

/3ao = p^K^T , (42) 

0ai = p^K^W , (43) 
A\ 2r n , 



/3 A2 = (— =-=zVX- (44) 

V V in J 3cpiV in 

with its essential nonlinear character manifest through its dependence on the average mag- 
netic induction B, the average current density J and the average Lorentz force J x B . 
One also observes that the Hall /3h requires a component of the current density fluctua- 
tions along the velocity fluctuations of the neutral fluid whereas the ambipolar effect thrives 
on the component of the current density fluctuations and the magnetic field fluctuations 
perpendicular to the velocity fluctuations in addition to the velocity- density correlation. 
We have used rigid or perfectly conducting boundary conditions (all surface contributions 
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vanish) while determining the averages. Here we consider what is known as the a 2 dynamo 

and take the mean flow Ve = 0. This actually determines the relative mean flow amongst 

the three fluids. The dynamo equation reduces to 
f)F} 

— = V x [aB - (3V x B] + r,V 2 B . (45) 
at 

The large scale magnetic field is written without the bar. The solutions of the dynamo 

equation for some representative cases were discussed in Krishan and Gangadhara ( 2008). 

Here the dynamo equation will be set up for the case where the neutral density fluctuations, 

through the ambipolar diffusion, are the major contributors to the a effect. Thus the dynamo 

equation becomes: 
f)Fi 

— = Vx [(a Aa - J x B)B - f3 A2 V x B]+ V W 2 B . (46) 

The equation even in this reduced form is still quite complicated. In order to have a glimpse 
of the role of the ambipolar term through the density fluctuations, we ignore the dissipation 
term and estimate the characteristic timescale of evolution of the field to be: 

t ~ 7 % • (47) 

where 



a M = ~ 2r -5 ^ P ' n ■ (48) 

The velocity-density correlation may be approximated by using the steady state mass 
conservation 

(7E + d,)(K + K')] = 0- (49) 



Taking the average and assuming \\V£ ■ Vp'J| ~ ||p^V • V^||, \\V n ■ Vp n || ~ ||PnV • V n \\ and 



||^-Vp^v n 'xVpUI; 

for the case of Ve = 0, retaining only the non-turbulent contribution, V n = — — 4^-. 
the characteristic timescale is found to be: 

t « rj^% . (50) 

T O 

cor rn 

where t A = w~ is the Alfven crossing time of the characteristic spatial scale L and V nA = 
is the Alfven speed of the neutral fluid. For the collisional plasma t A v in > 1 and t A is 
also expected to be greater than r cor as well. Thus depending on the ion to neutral density 
ratio the magnetic field evolution may take place on time scale as short as r cor . This rather 
crude estimate of the time scale is at best a pointer towards a possible rapid evolution of 
the magnetic field due to ambipolar diffusion on the solar atmosphere. 
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4 CONCLUSION 

In this first attempt at formulating the mean-field dynamo in a compressible partially ionized 
plasma, a complete set of multi-fluid- magnetic equations, in particular the Faraday law and 
the energy equations have been established with the hope that adequate computational 
resources would be deployed to understand this very important problem of the solar surface 
dynamo. Additional physics of the Hall and the ambipolar diffusion effects brings in a variety 
of spatial and time scales associated with the ion inertial scale, the gyrofrequencies and the 
collisional frequencies, the scales towards which the high resolution observations of the solar 
magnetic fields tend to lead us. 
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